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Recent experiments on fermions in synthetic gauge fields result in systems with a spin-orbit 
coupling along one spatial axis, a detuning field, and a Zeeman field. We show theoretically that the 
presence of all three results in interesting and unusual phenomena in such systems in the presence 
of a contact singlet attraction between the fermions (described by a scattering length). For two 
particles, bound states appear over certain range of the centre of mass momenta when a critical 
positive scattering length is attained, with the deepest bound state appearing at a nonzero centre of 
mass momentum. For the centre of mass momenta without a bound state, the gauge field induces 
a resonance like feature in the scattering continuum resulting in a large scattering phase shift. For 
many particles, we demonstrate that the system, in a parameter range, shows flow enhanced pairing, 
i. e., a more robust superfluid at finite centre of mass momentum. Yet another regime of parameters 
offers the opportunity to study strongly interacting normal states of spin-orbit coupled fermionic 
systems utilizing the resonance like feature induced by the synthetic gauge field. 
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Cold atomic systems with synthetic gauge fields [T] are 
expected to greatly enhance their ability to simulate ex- 
otic phases of matter from fractional quantum hall states 
to strongly coupled gauge theories. Spectacular exper- 
imental advances [2H1] have been successful in creating 
and studying systems with synthetic gauge fields. 

These advances have resulted in a flurry of theoreti- 
cal activity directed towards understanding interacting 
fermions in synthetic gauge fields. Uniform non Abelian 
gauge fields that produce a generalized Rasbha spin-orbit 
coupling have been investigated. Even with a weak at- 
tractive interaction, increasing the magnitude of such 
spin orbit coupling produces a crossover from a BCS state 
with large pairs to a BEC like state which is a conden- 
sate of new type of boson called the rashbon. [5] Super- 
fluidity of rashbons[H], BCS-BEC crossover with Zeeman 
fields etc[7H9], and transition temperatures of the rash- 
bon condensate [TOl |TT], have been studied (see [12], for a 
review) . 

Very recent experimental studies on fermions with syn- 
thetic gauge fields by the Shanxi[13j and MIT|14j groups 
have produced systems with a spin orbit coupling (A) in 
one spatial direction (equal mixture of Rashba and Dres- 
selhaus), a detuning (6) field, and a Zeeman (f2) field. In 
this paper, we explore the rich physics of such systems 
by theoretical studies of interacting fermions in them. 
We show that the simultaneous presence of the detun- 
ing and Zeeman terms (5 ^ 0, fl ^ 0), together with 
A, produces several interesting and unusual effects. At 
the two body level, a bound state exists only for a range 
of the centre of mass momenta, and the most strongly 
bound state (largest binding energy) appears at a finite 
centre of mass momentum. For centre of mass momenta 
where the bound state is absent, a resonance like feature 
with associated large change in the scattering phase shift 
appears in the scattering continuum. Interesting conse- 



quences that obtain in the many body setting include flow 
enhanced pairing (larger pairing amplitude for a flowing 
superfluid!). Our work also clearly bring out the pos- 
sibility of studying strongly interacting normal states of 
spin-orbit coupled fermions engendered by resonance like 
feature in the scattering continuum induced by the gauge 
field. 

The non-interacting part of the Hamiltonian of the sys- 
tem of interest is given by 



u 1 



(\k x - b)r* a , - Sir* 



(i) 

where cs are fermion operators, 5 and Vt are the de- 
tuning and two-photon Rabi coupling (Zeeman field) re- 
spectively, rs are Pauli matrices. As noted above, this 
type of spin-orbit coupling|15j has been realized in recent 
experiments. 0] [T3] E] The one-particle states of H have 
energies 



e a (k) = — - aB(k), B(k) = ^/ (\k x - 5)* + W , (2) 

where a = ±1 is the generalized helicity (simply called 
helicity below), with the associated states 
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where a (k) = arctan ( xk n _ & ^j + (1 — a)§. 

We introduce a contact singlet attraction between 
the spin- 4 fermions (such as that in a broad Feshbach 
resonance [16] ) 
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FIG. 1. Two body problem with no Zeeman field (6 7^ 

0,Q = 0): Singlet threshold e a th , two-body threshold Eth, 
bound state energy Ebs, and ground state energy e gs as a 
function of q. 
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FIG. 2. Singlet dos for 6 / 0, fi / 0: The singlet threshold 
el h coincides with the two-body threshold eth for all q. Plot 
shows g s (q,ui) for two values of q, q — A and q = — A. The 
vertical lines correspond to respective threshold energies. 
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where S+(g) = £ fe ^ ( 
and the bare parameter v is related to the s-wave scat- 
tering length a s through an ultraviolet cutoff A as = 
^ + A. We the study the ground state of a finite density 
of fermions in this system, after an investigation of the 
two-body problem. 

Two-Body Problem: It is convenient to introduce the 
following two body states: the singlet state 



and the state 



\q,k,s)=SHq,k)\0) 



\q,k,aP) =c\ +ka c\_ |0), 



(5) 
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where |0) is the fermion vacuum. We define the singlet 
amplitude as 



A a p(q,k) = (q,k,s\q,k,a/3), 
which for our system is 
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where 

Solution of ui of the secular equation 
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where V is the volume of the system, and 

£ a f}(q, k) = e Q (| + fe) + e /8 (- - fe), 

below the scattering threshold gives the energy of 
the bound state. The physics behind bound state 
formation [17] is determined by the singlet density of 
states (dos) 

g s (q,u) = ^Y. \A a p(q,k)\ 2 6{uj ~ e a p(q,k)). (11) 



The key aspect to be noted is that the singlet scatter- 
ing threshold, e% h (q), the smallest value of to such that 
g s (q,uj) = + , can be above the two-body threshold 
£th(q) = min fc e ++ (q, k), i. e., 



£th(q) < £ s th(q)- 



(12) 



We now discuss the solution of the two body problem 
in various regimes of parameters. We consider centre 
of mass momenta of the type q = qe x which contain 
interesting physics. 

Case A (5 ^ 0, ft = 0): Here we consider the case 
where A and 6 are non-vanishing. In the absence of the 
Zeeman field fl, we obtain an analytical expression for 
the threshold 



-th 
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with g s (q,u) = pJw^E^. The interesting aspect 
of this result is that, despite the presence of 5, which 
breaks the k x -H> — k x symmetry, the singlet dos satisfies 
g s (—q,uj) — g s (q,w). The singlet threshold coincides 
with the scattering threshold only in the regime q a < 
q < Qb (see Fig. [I]), where q a , qb depend on A and 6. For q 
outside this interval the scattering threshold Eth is below 
the singlet threshold. This is because the lowest two 
body state which corresponds to both the particles with 
+ helicity has no singlet component when q is outside 
the interval (q a ,Qb)- 

A scattering length a s > 0, results in a binding energy 
of l/al and a bound state energy, £& s (q) = —^2 + £ s th (q) 
(see Fig. [TJ. Note that the energy of this state is satisfies, 
£ bs{—q) = £bs(q), owing to the similar property of the 
singlet dos. This bound state is the ground state of the 
two-body system only in a regime of q, i. e., qi < q < 52 as 
shown in Fig.[T] where q\, q-i depend on A, 6 and a s . There 
is, therefore, a "first order" transition in the ground state 
at q\ (and 52) where the ground state abruptly changes 
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FIG. 3. Two body problem with Zeeman field (Q ^ 0): 

Singlet threshold ef h , two-body threshold e t h, bound state 
energy Ebs, and ground state energy e gs as a function of q. In 
this case the two-body threshold and singlet threshold always 
coincide. The deepest bound state (largest binding energy) 
appears at a nonzero value of q, at q m . Inset shows the q- 
asymmetry of the bound state energy. 



from the bound state with a large singlet component, to 
a "free" triplet state with both particles of + helicity. 

Case B (5 ^ 0, Vl ^ 0): How does a Zeeman field affect 
the picture just discussed? The key effect of the Zeeman 
field is on the singlet dos. In the presence of a Zeeman 
field however small, the singlet threshold e s th coincides 
with the two-body threshold e t h for all q. There is an im- 
portant additional effect. Although the singlet threshold 
is reduced, the singlet density of states does not have the 
characteristic Juj — e s th behaviour for all q. For the case 
shown in Fig. [2j for q = A, the singlet dos has the usual 
^Juj — e\ h threshold behaviour. In contrast, for q = —A 
the threshold behaviour is quite different. Analytical cal- 
culations show that the singlet dos in this case goes as 
(cj — e^J 3 / 2 . It can in fact be shown that this is the 
case whenever q is outside an interval (q a ,qb) (Fig. [3]). 
Indeed in Fig. [2] q = — A corresponds to a q that lies 
outside the said interval. The gist is that when 5^0 
and f2 7^ 0, the density of states is not symmetric in q, 
i. e., g s (—qe x ,ui) ^ g s (qe x ,ui). This is to be contrasted 
with the previous discussion when Q = 0, where g s was 
g-symmetric even when 5^0. 

These characteristics have interesting effects on the 
ground state of two particles. A positive critical scat- 
tering length is necessary to produce a bound state, ow- 
ing to the Zeeman field. Even for a positive scatter- 
ing length smaller than this critical scattering length, as 
shown in Fig. |3j a well defined bound state exists only 
when qi < q < qi- The deepest bound state appears at a 
finite momentum q m as shown in Fig. [3j and the bound 
state energy is not at symmetric function of q (Fig. [3jin- 
set). 

When q is outside the interval q\ < q < (& (Fig. [3]), 
bound states cease to exist. However, the "remanent" of 
the bound state appears as a "resonance like" feature in 
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FIG. 4. Scattering phase shift: For q ^ (91,92 ) (see Fig. [3), 
a resonance feature appears in the scattering continuum with 
a concomitant large increase in the scattering phase shift 77. 



the scattering continuum with a large change in phase 
shift 77 |18| appearing at a g-dependent energy in the 
scattering continuum (see Fig. B. This is akin to what 
is found in narrow Feshbach resonances (without spin- 
orbit coupling) 19J. The notable feature here is that this 
resonance like feature is induced by synthetic gauge field, 
i. e., the spin-orbit coupling along with the detuning and 
Zeeman fields, even when the interaction between the 
fermions is from a wide resonance described by an energy 
independent scattering length. As is evident, this will 
have profound consequences for the many body problem 
as strong fermion-fermion interactions can be produced 
in the scattering continuum. 

Many-Body System: We now consider a system of 

it 3 

^ with an associated 



fermions at a finite density po 

energy scale Ep = . As was learnt from the analysis of 
the two body problem, interesting physics occurs when 
both 6 and fi, along with A, are non-zero, and we shall 
only consider such a case here. 

Noting that the bound state is deeper at a finite q (see 
Fig. [3]) , we consider many-body states corresponding to 
a "flowing superfiuid" . This can be done by constructing 
a mean field theory with v(S^(q)) — After some 

algebra, we obtain the following mean field Hamiltonian 



Umf = ^ t (g,fc)e(q,fc)*( g ,fc) 



5>(ff(^-fc)-M 
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: 

V 



(14) 



where the symbol stands summation over all k with 
k x > 0, *(fc) 



( c (|+fe)t c (i+k)i c (f-fc)T c (f-fcu)' 



W(q,k) 



H(Z +k)- fj.1 
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FIG. 5. Dependence of pairing amplitude on q: Landau 
critical momentum (q at which A vanishes) for flow + direc- 
tion is larger than in — direction. Inset shows flow enhanced 
pairing where A attains a maximum value at a nonzero q. 



where H(k) is the matrix defined by eqn. |l]), 1 and cr y 
are the identity and Pauli-y matrix, all three of which 
are 2x2 matrices, and fj, is the chemical potential. The 
thermodynamic potential P at a temperature T is 



VP(T, n) = -TJ2 H 1 + e- E ^ q ^/ T ) 



q ,\ VA 2 



(16) 



$>(ff(|-fc)-/xi 



ka 



where n runs over the four eigenvalues of H(q, fe). Equi- 
librium values of A and [i are obtained by the gap and 
number equations derived from eqn. (16). 

Let us discuss the results of the calculations. The first 
point to be noted is that in the presence of the Zeeman 
field, any attractive interaction does not produce a su- 
perfluid state. A non-vanishing critical scattering length 
a sc (which is q dependent) is necessary to induce pairing. 
This is analogous to what is known in the usual Fermi 
gases (without spin-orbit coupling) with imbalance. The 
physics of this, again, owes to the singlet density of states. 

Fig. [5] shows the dependence of A on q (q = qe x ). 
The magnitude of q at which A vanishes (related to the 
Landau critical velocity) is different for flow in the + 
direction and — direction. In fact, the critical velocity 
for flow in the + direction is about 5% larger than for 
flow in the — direction. This asymmetric dependence of 
A on q is already seen near q — (see inset of Fig. [5]) , 
where it is clearly seen that A peaks at a non-zero value 
of q which is an instance of flow enhanced pairing. We 
emphasize again that this phenomenon is absent when 
5 ^ and Cl = 0. 

The results just discussed may produce an impression 
that the asymmetric nature of the critical velocity on the 
direction of flow and the flow enhanced pairing, albeit 
clearly demonstrated, is a small effect that might be dif- 
ficult to detect in experiments. While this is true at large 
densities, the effect is strongly enhanced at lower densi- 



ties as we show by the analysis below. Consider a system 
where 6/A 2 and fi/A 2 are kept fixed, and kp/X <C 1, i. e., 
a low density system. Analysis of q dependent pairing 
in this case is simplified by transforming the mean-field 



Hamiltonian (eqn. (14)) to the helicity basis 
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(17) 



where £ a (k) — e a {k) — fj,. Setting T = 0, we focus on 
the regime of q in (91, 92) (see Fig. [3]). For a scattering 
length such that — > — , in the limit of hp/ A -C 1, we 
have |A| -C \fj,\. Analysing the gap equation in this limit 
via a second-order perturbation theory in A gives 



ka/3 



\A a0 {q,k)\< 



1 ~ 



\A a p(qM 2 
(18) 



which is the secular equation (eqn. Q) corresponding to 
the two body problem. What we have demonstrated is 
that in the regime kp /A -C 1, and q € (51,92), the sys- 
tem is a BEC of the tightly bound pairs described by the 
bound state shown in Fig. [3] What is even more inter- 
esting is that, this BEC becomes a "more robust" super- 
fluid with increasing q starting from q\ ; indeed, there is a 
9m € [q\, (fe) where the binding is the largest (see Fig.|3|, 
clearly demonstrating that flow promotes superfluid pair- 
ing in this regime. It will be interesting to explore ways 
to demonstrate this in experiments. 

The analysis above bring out another very interesting 
aspect of this system. What happens in the regime where 
hp /A <C 1 and q ^ ((71,92)? The system will be a normal 
fluid. What is remarkable is that this normal fluid will 
be strongly interacting due to the presence of the reso- 
nance like feature in the scattering continuum. This state 
should be accessible in the current experiments. The idea 
would be to tune the parameters to a regime where qi > 
so that even at q = the system is in this strongly inter- 
acting normal state with the Fermi energy tuned to the 
position of the resonance to produce strong interactions. 
In fact, the broad resonance of Li at 830 Gauss can be 
utilized to study such strongly interacting Fermi liquids, 
which arc likely to show strong pseudogap features. An 
all important question is weather such a state is unstable 
towards other orders such as magnetism. Our findings, 
therefore, provide an interesting direction of future re- 
search. 

In summary, we have demonstrated that spin-orbit 
coupled systems, with detuning and Zeeman fields, pro- 
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vide scope for study of a variety of new phenomena of 
interest to condensed matter physics. We hope that this 
work stimulates further experiments. 
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